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Zigzag line defects and manipulation of colloids in 
a nematic liquid crystal in microwrinkle grooves 

Takuya Ohzono 1 & Jun-ichi Fukuda 1 



Spatially confined liquid crystals exhibit non-uniform alignment, often accompanied by self- 
organised topological defects of non-trivial shape in response to imposed boundary conditions 
and geometry. Here we show that a nematic liquid crystal, when confined in a sinusoidal 
microwrinkle groove, exhibits a new periodic arrangement of twist deformations and a zigzag 
line defect. This periodic ordering results from the inherent liquid crystal elastic anisotropy and 
the antagonistic boundary conditions at the flat liquid crystal-air and the curved liquid crystal- 
groove interfaces. The periodic structure can be tuned by controlling the groove geometry 
and the molecular chirality, which demonstrates the importance of boundary conditions and 
introduced asymmetry for the engineering of topological defects. Moreover, the kinks in the 
zigzag defects can trap small particles, which may afford a new method for manipulation 
of colloids. Our system, which uses easily fabricated microwrinkle grooves, provides a new 
microfabrication method based on the arrangement of controllable defects. 
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Nematic (N) liquid crystals (LCs) are anisotropic fluids with 
long-range molecular orientational order and no positional 
order of their constituent molecules 1 . The local orienta- 
tional order of an NLC is specified by a unit vector, n, known as the 
director. The free energy of an NLC in terms of n, called the 
Frank elastic energy is expressed by F=(l/2)Jdr[K n (divw(r)) 2 + 
K 2 2(n(r)-cm\n(r)) 2 + K 33 (n(r)xcm\n(r)) 2 ], where K n , K 22 and K 33 
are elastic constants for splay twist and bend distortions, respec- 
tively. The static profile of the director is determined by minimising 
the Frank elastic energy under the boundary conditions imposed 
at the surfaces confining the NLC. The geometry of the confining 
surfaces and the interaction between a surface and NLC molecules, 
called the surface anchoring, have a crucial role in the control of LC 
alignment. Although a uniform n would minimize the Frank elas- 
tic energy density in the bulk, various distorted alignments appear 
when the LC is confined in specific geometries, such as the twisted 
N cell, hybrid alignments 2 " 7 , capillary tubes 8-10 , emulsion drop- 
lets 11 ' 12 , thin cells 13 and shells 14 . Patterned surfaces 15-17 are also 
used to induce regularly distorted alignments. The resulting distor- 
tions, which are often periodic or include topological defects, make 
LCs useful for electro-optical applications, assembling colloids 18 ' 19 
and for the study of topological problems in a vector field 20,21 . 

In a typical hybrid-alignment system, where one of the two con- 
fining surfaces imposes normal anchoring, and the other imposes 
planar anchoring, a topological defect line, known as a disclination, 
generally forms during the equilibration and is stabilized by a non- 
flat boundary 7 . When the director n is perpendicular to a nearby 
disclination line, the disclination is called a wedge disclination, and 
the winding number is well defined; in this case, it is ± 1/2. A half- 
integer winding number is allowed by the equivalence of n and -n 
in an NLC. 

In this paper, we consider an NLC in a microgroove, where the 
NLC aligns normal to the flat air-LC interface, and the curved LC- 
groove interface imposes a planar alignment with its easy axis per- 
pendicular to the groove direction. We utilize a recently developed 
method to manipulate liquids into microfilaments in the curved 
grooves of self-organized micro wrinkles 22 . We have discovered a 
novel symmetry-broken, self- organized structure with periodic 
domains of alternating twist directions and a zigzag disclination line 
with a winding number of + 1/2. The geometrical arrangement in a 
previous study reporting a zigzag disclination 9 ' 10 is entirely differ- 
ent from ours; a flat capillary tube employed there does not allow 
the NLC to come into contact with air, and the surface of the tube 
imposes normal alignment. Moreover, our microwrinkle system 
allows the geometrical parameters, such as the groove wavelength 
in the micrometre range, to be easily controlled. We propose that 
the non-trivial configuration arises from the intrinsic anisotropy 
in the elasticity of an NLC with K n ^K 22 ^K 33 , and that a smaller 
K 22 yields a spontaneous twist deformation below the disclination, 
which rotates it away from the groove direction. 

Results 

Observation of twisted domains with zigzag line defects. Our 

microwrinkle surface 22-25 (see Methods) has a sinusoidal undulating 
surface with a specific wavelength A, and a depth D; the grooves can 
serve as open capillaries for a liquid 22 . The liquid fills about three- 
quarters of X (Fig. la, Supplementary Fig. SI) and is called a liquid 
filament. A typical NLC, pentylcyanobiphenyl (5CB), was confined 
as liquid filaments. The air-NLC interface is almost flat in contrast 
to the bottom- curved surface of the wrinkle groove. The maximum 
thickness of the liquid filament at the centre, h, is slightly smaller 
than D. 

The air-NLC interface imposes normal alignment, where the tilt 
angle measured from the z axis 6~ 0, and the azimuthal angle meas- 
ured from the x axis, 0, can take any value (Fig. la). In contrast, the 
wrinkle surface imposes planar alignment, where 6~7i/2 and 0-0, 



or 7T, with the easy axis direction perpendicular to the groove direc- 
tion (Supplementary Fig. S2). The local easy axis is in the x-z plane 
and along the tangential direction of the sinusoidal curve. 

Symmetry-broken, periodic optical rotations were observed 
along the groove direction with a periodicity A (Fig. lb-e). Fig. If 
shows the proposed alignment structure from the analysis of the 
optical images (Supplementary Fig. S3); 6 decreases from nil to 0 
from the bottom to the top of the groove; therefore, the optical rota- 
tion, which reaches a maximum in the centre of the domain, can be 
attributed to twist distortions in the bottom half, 0~ ±25°. The black 
zigzag line that deviates from they axis in Fig. lb,d is a disclination 
line that changes its direction at the domain boundaries, separat- 
ing domains with opposite twist handedness. A different observa- 
tion of an inclined sample shown in Supplementary Fig. S4 clearly 
demonstrates that the winding number of the disclination line is 
positive, + 1/2. Obviously, with our concaved wrinkle groove with 
planar anchoring, a + 1/2 disclination line involves smaller splay- 
bend distortions than a - 1/2 one. The angle of the disclination line 
from the y axis, Q~ ± 8° and (j)Q > 0, suggesting that the local direc- 
tion of the disclination line and the twisted local director tend to 
make a smaller angle rather than retain a right angle (the sign of Q 
is defined so that a clockwise rotation gives a positive value of Q; see 
Fig. If). This periodic domain structure is stable at room tempera- 
ture (-300 ± 5 K) for at least a month, unless the LC evaporates. 

The periodic domain structure is not observed when X is below 
-1.5 Jim (Fig. lg), or /z<150nm. If h is much smaller than the 
characteristic length associated with surface anchoring, R C ~K/W 
~l|im (refs 11,21) (iO10 -11 N, and the anchoring energy 
W~10 -5 Nm -1 ), the Frank elasticity dominates the surface 
anchoring, which results in a more uniform NLC alignment without 
a + 1/2 disclination. The periodic domain structure is highly associ- 
ated with the presence of a disclination line; thus, the formation of a 
periodic domain structure is suppressed when h or X are small. 

Breaking symmetry of periodic domain structures. Next, we dis- 
cuss several ways of breaking the symmetry between two neighbour- 
ing domains with opposite twist handedness and the consequences 
of introducing asymmetry. The twist handedness can be determined 
by introducing a chiral reagent, such as CB15, that induces right- 
handed twist (Fig. 2a-c). The liquid filament of the chiral nematic 
(N*) with a right-handed twist distortion does not form a periodic 
domain structure (Fig. 2b, bottom). Immediately after the fusion of 
the N and N* liquid filaments, where the CB15 density gradient still 
exists, the fused liquid filaments show a continuous transition from 
the periodic domain structure to the right-handed twist region with 
no specific discontinuities (Fig. 2c). This demonstrates that N* sta- 
bilizes one of the twisted domains in the periodic domain structure. 
The symmetry breaking also occurs when the angle between the 
easy axis of the NLC alignment and the groove direction deviates 
from a right angle, which is often observed at microwrinkle dis- 
locations (Fig. 2d-e). In this case, one of the rotated disclination 
line segments is stabilized over the other because of the symmetry 
breaking, and thus, the formation of a periodic domain structure is 
suppressed. In both cases, one prescribed twist handedness domi- 
nates, indicating that the periodic domain structure can be con- 
trolled by changing the chiral strength of the NLC or the patterned 
alignment conditions. 

Manipulation of colloids using zigzag defect lines. Disclination 
lines can trap small particles to reduce the total elastic energy of the 
system 26-28 ; the particles remove the energetically costly disclina- 
tion cores. We indeed observed silica beads trapped at the disclina- 
tion line at which the strongest elastic deformation is likely (Fig. 3). 
It should be noted that the particles cannot move freely along the 
disclination line, but are trapped selectively at its kinks. The kinks in 
the disclination line are presumably more energetically unfavourable 
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Figure 1 1 Periodic alignment with a zigzag line defect in a NLC confined in microwrinkle grooves, (a) Geometrical confinement and hybrid anchoring 
conditions for a liquid filament of an NLC on a microwrinkle groove, where D is the maximum depth of the groove, h is the maximum thickness of the 
liquid filament, X is the wavelength of the microwrinkles, and w is the width of the liquid filament. The blue curves indicate the profile of the sinusoidal 
groove. The definition of the polar and the azimuthal angles 6 and <j> of the local nematic director n (bold blue line) is also given, (b-e) Optical images with 
different configurations of the polariser (P) and analyser (A). Scale bar is 20um In (e), a sensitive tint plate is used with the configuration shown as S. The 
colour changes depending on the angle between the local n and the direction of S; blue (orange) colour suggests that n is parallel (normal) to the direction 
of S. The nail symbols (T) indicate the projection of n onto the plane of the page (the head of the nail comes out of the page), (f) Schematic of the director 
structure n(r) with the definitions of the maximum twist angle <j> and the rotation angle of the disclination line (red line) from the y axis Q. A is defined as 
the wavelength of the periodic domain structure. The red dots indicate disclination kinks. The blue (orange) colour indicate domains with (f)< 0 (</>> 0), 
which corresponds to the colours in (e). The domain boundaries correspond to the black inclined line passing through the centre line of the groove in (c). 
(g) Plot of A versus X; mean + s.e.m. are shown and the number of samples for each point is over 50. The data at X = 2\im show larger error relative to the 
mean value, because the periodic domain structure is less regular. A periodic domain structure was not observed with X below -1.5 Jim. 



than the straight parts, which result in stronger binding of beads at 
the kinks. In addition, when the system is heated to the isotropic 
phase, the beads are no longer trapped and show random motion 
(Supplementary Movie 1). Our observation demonstrates that top- 
ological defects can have an active role in the control of colloidal 
assembly in LCs, whereas in usual LC colloids, topological defects 
are formed passively in response to colloid- surface anchoring 19 . 
Zigzag disclination lines can thus provide a novel way of manipulat- 
ing colloidal particles. 

Discussion 

We discuss why a disclination line in a narrow groove can take a 
direction oblique to the groove, resulting in a zigzag shape. Fol- 
lowing the discussion by Buscaglia et at. 29 on the spontaneous 
twist deformations in an NLC under a straight disclination line, we 
assume that the dominant contributions to the total elastic energy 
of an NLC in the groove are from the elastic deformation in the 
cylindrical region around the disclination line, and from the twist 
deformation in the NLC below the cylindrical region. In a simpli- 
fied case where K\ i = K33 = K, the total elastic energy per unit groove 
length is given by: 



a c + J cos 2 (Q + 0) + ^ sin 2 {Q + 6) 1 , 

^(ftfl-K-E-V - ^ - ^70 2 - (1) 

cos L2 2 

The first term is the contribution from the cylindrical region around 
the disclination line, where K^(K/4)K\n(R/r c ), R is the radius of the 
cylinder, r c is the microscopic cut-off, and 0 is the azimuthal angle 
of n at the bottom of the cylindrical region. The parameter a c , which 
is of the order of 0.1, arises from the contribution inside the micro- 
scopic cut-off where the continuum Frank elasticity no longer holds. 
The second term accounts for the twist deformation below the cylin- 
drical region. The symbol /is a monotonically increasing function 
of the anchoring strength at the groove surface, and takes a value 
between 0 and K22 A/d 2 , where A is the area of the region where the 
twist deformation is relevant, and d is the distance between the bot- 
tom of the cylinder and the groove surface. Equation (1) indicates 
that a positive (j) yields a smaller £ to tal for 12 > 0 and vice versa. This is 
in agreement with our experimental findings described above. See 
Fig. 4a for the geometry, and Methods for theoretical details of the 
above and the following discussions. 
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Figure 2 | Breaking the symmetry between two neighbouring domains with opposite twist handedness, (a) Schematic of the fusion process of liquid 
filaments protruding from the NLC and NTLC droplets. The thick arrows indicate the direction of the capillary action, (b) Optical image of the liquid 
filaments before fusion. Note that the liquid filaments of N* do not form a periodic domain structure (bottom of (b)). Scale bar is 30jim. (c) The fused 
liquid filaments of N and N* with the gradient of the chiral reagent. The nail symbols (T) indicate the projection of n onto the plane of the page (the head 
of the nail comes out of the page), (d) Optical image of the periodic domain structure in a large area that includes randomly produced dislocations of 
microwrinkle grooves. Scale bar is 50 jam. (e) Magnified image of the region indicated by the dotted rectangle in (d) with a slightly rotated polariser (P) 
and analyser (A) to differentiate the domains with different twist handedness. The groove direction deviates slightly from the y axis around the dislocation, 
as indicated by the yellow arrows. The absence of periodic domain structures is highlighted by the white arrows. The easy axis is in the horizontal 
direction. Liquid filaments are not formed at the dislocation, because the shallow grooves destabilize the liquid filament 22 (dotted ellipse). 
Scale bar is 20um 



a bed 




Figure 3 | Trapped silica beads at the domain boundaries of a periodic domain structure, (a) Schematic of the trapped beads at the kinks of the zigzag 
disclination line, or domain boundaries of the periodic domain structure. (b,c) Optical images of the periodic domain structure with two trapped silica 
beads at the disclination line kinks. The arrows indicate the positions of the beads. In (b), a tint plate (S) was used. The colour changes depending on the 
angle between the local n and the direction of S; blue (orange) colour suggests that n is parallel (normal) to the direction of S. (d) Optical image of the 
released silica beads just after the transition from a nematic to an isotropic phase (see also Supplementary Movie 1). The disclination line disappears and 
the beads are free to move in the isotropic liquid filament. The same part of the system is pictured in (b,c and d). Scale bar is 15jim. 



A straight disclination line along the groove (Q = 0) becomes 
unstable when: 



c)(l-») 



(2) 



In most of the NLCs of rod-like molecules, K 2 2 < Kn < K33 (ref. 30) 
and thus the right-hand side of equation (2) is positive. The numeri- 



cal analysis of equation (1) for various values of K22/K and y/K, 
given in Fig. 4b -d, shows that Q that gives the minimum of £ t otal 
(referred to as Q m later) is largest at an intermediate value of y/K, or 
an intermediate anchoring strength at the groove surface, and that 
Q m can be larger for smaller values of K22/K. 

Note that a disclination line with Q = 0 cannot be unstable when 
K 2 2=K (so-called one-constant elasticity), or 0 is rigidly fixed to zero 
(in both cases the first term of eq. (1) is a monotonically increasing 
function of Q at 0 < Q < nil). One can understand that this instabil- 
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Figure 4 | Geometry of the system for the theoretical arguments and 
numerically evaluated free energy of the system, (a) Schematic of the 
geometry of the system. R is the radius of the cylindrical region involving 
a disclination line. r c is the cut-off radius. The director profile n is given by 
nails in red (the head of the nail comes out of the page). When the twist 
angle (p=0,n lies in the x-z plane, or the plane of the page (left), whereas 
non-zero (p yields n out of the x-z plane (right). In the region shaded in dark 
yellow, with area A and height d, the twist distortion leading to equation 
(8) is present. The region shadowed by slashes is the bottom of the groove, 
(b-d) Rescaled energy of the system E {Q{a \(£2)/K as a function of X2for 
K 22 /K= 0.25 (b), 0.5 (c), and 0.75 (d). Here, a c in equation (1) is set to 0.1. 
The numbers are the rescaled anchoring strength y/K. The points indicate 
the minimum of E t0{a \(£2)/K. The numbers in bold are the smallest y/K 
violating equation (2). 



ity occurs because the elastic energy of the disclination line, the first 
term of equation (1), is reduced by adopting positive 0, with the 
expense of the twist deformation below it (the second term of equa- 
tion (1)). For this mechanism to be possible, the twist deformation 
cannot be energetically costly, requiring relatively small K 22 and 
thus small y/K. However, when yis so small (typically in the case of 
weak surface anchoring) that the second term of equation (1) is not 
effective, to minimize the first term, 0 tends to be close to nil - Q. 
Then the denominator of the first term dominates the behaviour 
and again Q m ~ 0, as seen from Fig. 4; the instability is effectively 
invisible. Thus, the instability of a disclination line along the groove 
direction resulting in finite Q is a consequence of a delicate reduc- 
tion of the disclination energy caused by finite (f). 

The disclination line takes alternate directions Q m and -O m 
along the groove, both of which are allowed by symmetry, so that the 
disclination line can avoid strong elastic distortions of the director 
when it approaches the groove walls. This alternation results in the 
experimentally observed zigzag form. The wavelength A of the zig- 
zag can be roughly estimated by using the elastic energy associated 
with the approach of disclination line to the groove walls, and the 
elastic energy in narrow regions of width from the abrupt change 
in direction of the disclination line. The former free energy, per unit 
length of the groove, is estimated as/i — K(h/X)(A tan£2 m M) 2 , and 
the latter energy from the change in direction is roughly estimated 
as/ 2 — Klrf/A. Numerical prefactors of order unity have been omit- 
ted. When/i + f 2 is minimized with respect to A, Aoc(tanI2 m )~ 2/3 
is obtained, and A<^l(l^/h) l/3 . The aspect ratio of the groove k/h 
is fixed experimentally to be almost constant 22 . If it is assumed 
that Itf^A, as expected in an NLC under strong surface anchoring 
where only the characteristic lengths of the geometrical confine- 
ment (X and h) determine the characteristic length scale of the 
elastic distortion (/j), then A<^X, in agreement with experimental 
observations (Fig. lg). 

The self-organized topological structure of an NLC with a zigzag 
disclination line reported here is an interesting example of spon- 
taneous symmetry-breaking phenomena arising from the coupling 
between material anisotropy, and asymmetric boundary condi- 
tions in a confined geometry. It may provide a new experimental 
system for studying basic and unexplored properties of topological 
defects not only in various LC phases, but also in other condensed- 
matter systems with a vectorial order. Our experiments show that 
our easy-to-make self- organized microwrinkles are quite useful 
for engineering the new zigzag disclination line structures and for 
manipulating small particles. The tunability of the shape of micro- 
wrinkles by simply exerting strain 22 ' 25 allows further dynamical 
control of the defect structure. Other efficient ways of controlling 
the defect structure would be the application of an electric or a mag- 
netic field, and the tuning of LC elasticity by varying temperature 
or by changing the chirality and twisting power of the LC material. 
The variety of such technical advantages facilitates the microman- 
ipulation of LC defect structures at a small scale with a very simple 
procedure. Our results demonstrate a new possibility of fine control 
of colloidal assembly without relying on conventional top-down 
lithographic techniques, and further knowledge of how the surface 
anchoring and the size of colloidal particles affect the trapping phe- 
nomena will sophisticate the manipulation of colloidal assembly 
using self-organized defects. Thus, our system with controllable 
geometry provides a route to a novel microfabrication method 
for colloidal assembly, optical patterning and guided polymeri- 
zation and crystallization 31 , with its self- organized structure as a 
template. 



Methods 

Preparation of microwrinkles. Microwrinkles with aligned grooves were 
prepared as follows, which is identical to the method reported by Ohzono et al. 22 
A polydimethylsiloxane (PDMS) elastomer (Sylgard 184, Dow-corning) was cured 
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on a smooth Si wafer. The smooth surface of the PDMS substrate (12x12x6 mm 3 ) 
was weakly treated by an Ar plasma (SEDE-P, Meiwa Forsis) to enhance the wet- 
ting of N-methylpyroridone (Sigma- Aldrich). Soon after an N-methylpyroridone 
solution of poly-amic acid (Poly(pyromellitic dianhydride-co-4,4'-oxydianiline), 
Sigma- Aldrich) was coated on the activated PDMS surface using a spin-coater at 
5000r.p.m. for 90s (MSA- 150, Mikasa), the substrate was uniaxially compressed 
with the strain of approximately - 0.04 ( - 4%). The sample was heated under 
compression at 180°C for 3 h to induce the transformation of poly(amic acid) to 
polyimide. During the temperature increase in this heating process, the PDMS 
sample dilated to the direction perpendicular to the compression due to thermal 
expansion. Thus, the main chains of the polymer plausibly aligned anisotropically, 
which imposed the observed alignment of the cyanobiphenyl NLCs. After cooling 
down to 25 °C and the release of the compression, aligned microwrinkle grooves 
were formed with the easy axis of the LC alignment perpendicular to the groove 
direction. The wavelength of the groove A, 0.6-20 \im, was controlled by varying 
the concentration of the poly-amic acid solution from approximately 1 to 10 wt%. 
The aspect ratio of the microwrinkles, D/A, where D is the groove depth, was 
almost constant, 0.13 ±0.01. 

Formation of liquid filaments of LCs. Liquid filaments were formed on the 
microgrooves, by placing a droplet of an NLC with a sub -millimetre diameter on 
the wrinkled surface. With the present aspect ratio, 5CB was spontaneously sucked 
into the microgrooves via capillary action 22 . Alternatively, the liquid filaments were 
formed by simply coating microwrinkles with 5CB (see Supplementary Fig. SI 
for further details). The LC alignment was observed using a transmitted polarized 
microscope (BX-51P, Olympus). To differentiate small optical rotations visually by 
the colour difference, a 530-nm sensitive tint plate (1 A, U-TP530, Olympus) was 
used as a test plate compensator under crossed-Nicol conditions. A N* phase doped 
with a right-handed chiral reagent, 4-cyano-4'-((S)-2-methylbutyl)biphenyl (CB15, 
Merck, ~3 wt%) was used to investigate the effect of chirality. Droplets of N*and N 
(pure 5CB) LCs were placed a certain distance apart (-100 Jim), and the liquid fila- 
ments protruding from each droplet made contact on each microgroove and fused 
into liquid filaments connecting the droplets. A dispersion of melamine- coated 
silica beads with an average diameter of -600 nm (Optbeads, Nissan Chemical, 
-0.5 wt%) in 5CB was used for trapping colloids. The maximum thickness of the 
liquid filament, h ( > 1.5 urn), was large enough for the silica beads to be fully im- 
mersed in 5CB without touching the bottom groove surface. We used a transmit- 
ted polarized optical microscope to observe the beads trapped at the domain 
boundaries. 

Details of the theoretical arguments. For the derivation of equation (1), let us 
first explain the free energy of the disclination line evaluated in a cylindrical re- 
gion of radius R. See Fig. 4a for the geometry of the system. We consider a straight 
disclination line uniform along the y direction. When the disclination is a wedge 
disclination with the director n lying on a plane perpendicular to the disclination 
line (an x-z plane), n = (sin 6, 0, cos 0), where 0 is a function of x and z. Now we 
are interested in the effect of the twist deformation on the energy of a disclina- 
tion line, which we assume can be incorporated by a single variable (p describing 
the rotation of n along the z axis. We further assume that (p is independent of 
x, y and z. Then, the director field can be given by n = (sin 6 cos (p, sin 6 sin (p, 
cos 6). The Frank free-energy density /prank> under the simplifying assumption of 
K± i = = K, is 



/Frank ~ ~ K 



(3) 



with x = xl yj cos 2 (p + (K 22 I K) sin 2 (p. A familiar solution involving a disclina- 
tion line of strength +1/2 that minimizes the Frank elastic energy Jdr/p ran k is 
6 = 00 + (1 / 2) arctan(z / x) with 0 0 being an integration constant, and thus 
/Frank = (K / S)(x 2 + z 2 )~ l . The free energy of a disclination line per unit length 
is then expressed as 



^disclination = E c + ~ J dr j Q d X 



^cos 2 (p + (K 22 1 K) sin 2 q> 



+ sin x 



- E c + f In — j a/cos 2 (p + (K 22 I K) sin 2 <p 



(4) 



where we have introduced a transformation of variables x = rcos% and z = rs'm%, 
together with a microscopic cut-off r c . The contribution E c , of the order of K 
(ref. 1), comes from the region r < r c in which the description of the orientational 
order in terms of n is no longer valid. In the case of our interest, \n(R/r c ) — 10 
(ref. 1), and thus E c = a c K"with a c — 0.1. When the disclination is not oriented along 
the y axis (groove direction) and the angle between the disclination and the y axis 
is Q (Fig. If), the free energy of the disclination line per unit length along the 



groove direction is £disclination/ cos ^- Moreover, when the azimuthal angle of n, the 
angle between n and the x axis, is (f) at the bottom, (p = Q+(f). Thus, 



^disclination . 

cos 12 



a c + ^cos 2 (Q + (/)) + K ^ sin 2 {Q +"<; 
cosX2 



(5) 



with K={nKlA)\n(Rlr c ), yielding the first term of £ t otal> equation (1). We note here 
that the form of the disclination energy presented in ref. 29 is different from ours, 
equation (4). This is because Buscaglia et al. 29 made an oversimplifying assump- 
tion for the profile of n around a disclination; in our notation, they employed 
6= ( l/2)arctan (z/x), valid only when K x 1 = K 22 - ^33, and applied the profile to the 
general cases of Kn^K 22 ^K^. 

In the evaluation of the elastic energy due to the twist distortion below the 
cylindrical region containing the disclination line, we make a crude assumption 
that the twist distortion there is present in an area of A at a cross section parallel to 
the x-z plane, and that / is the width of the groove bottom in contact with the NLC. 
The azimuthal angle of n at the groove bottom is referred to as <py. Then the free 
energy associated with the twist distortion and that with surface anchoring at the 
groove bottom are, per unit length along the groove direction, 



B twist — ~ K 22 A 



1 2 
^anchoring -~Wlsm $ 



(6) 



(7) 



In equation (6), d is the distance between the bottom of the cylindrical region con- 
taining the disclination and the groove bottom, throughout which the twist distor- 
tion is present. In equation (7), W implies the anchoring strength, and it is assumed 
that I (j)y I <C 1. By substituting the solution of d(-Etwist + -^anchoring) / ^^b - 0 as an 
equation for tyb into equations (6) and (7), one can eliminate (py, and 
obtain 



^twist + -^anchoring — ■> 

constituting the second term of £ t otal> equation (1), with ^(K 2 2WlA/d 2 )/(K 2 2A/ 
d 2 + Wl). One can easily see that 7 is a monotonically increasing function from 
0 to K 22 Ald 2 with increasing W from zero to infinity. 

For the derivation of equation (2), we expand equation (1) up to second 
order in Q in the numerator of the first term, and in Q + (f) in its denominator. 
The result is 



c(2) . 
c total ' 



-(l + a c )Q 2 +-(^- 
2 2V K 



1 + < 



1 ^2 



(9) 



Here the superscript (2) is to emphasize that the expansion is kept up to the second 
order. By eliminating (j) after solving 5-E^otal / - 0 as an equation for </>, and 
substituting the solution to equation (9), one obtains 



c(2) . 
D total ' 



l + a r 



K 22 
K 



y-K(l 



K22 , 

K > 



(10) 



Thus, it can be shown that Q = 0 becomes unstable when equation (2) is fulfilled. 

To find Q that minimizes £ t otal (referred to as Q m later), numerical calculations 
are necessary. In Fig. 4b-d, we plot E tota \(Q) for several values of ylK and K 22 /K 
after the elimination of 0 by solving 3-Etotal / ^0 = 0 and substituting the solution 
to equation (1), the same procedure as in the above discussion. One can find that 
the range of y/K"in which Q m is non-zero becomes smaller with the increase of 
K 22 /K towards 1 (as also seen from eq. (2)). It is also clearly seen that when ylK 
approaches zero, corresponding to the limit of weak anchoring, Q m approaches 
zero as well, and that the largest Q m for fixed K 22 /K, which becomes larger with the 
decrease of K 22 /K, is attained at some intermediate value of y/K", or intermediate 
surface anchoring at the groove bottom. 

Finally, we estimate j\ necessary for the evaluation of A. The elastic energy 
associated with the approach of the disclination line to the groove walls, per unit 
length of the groove, is estimated as 



fx{y)^-c x Kh\ 1 + 1 

J1/ 2 1 \Xl2-u{y) X/2 + u(y) 



(11) 



where u{y) is the displacement of the disclination line from the centre of the 
groove, and q is a numerical constant of order unity. Equation (1 1) is composed 
of the elastic energies in a region between the disclination line and one of the 
walls, involving the distortion of the director n with a characteristic length 
A/2-w(v), and in another region between the disclination line and the other wall 
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whose distance is A/2 + u(y). The assumption that the disclination line is made 
up of a straight line segment u{y) =y tan Q m at - A/4 <y < A/4 yields the spatial 
average of/i as 

A/2J-A/4 /n 

2c,Kh , 2A + AtanX2 M 

= l - In— — (12) 

AtanX^ 2A-AtanX^ v ; 

3A I 2/1 J ' 

resulting in the form of/i given in Discussion apart from a numerical factor 
(in Discussion, the overbar of f± has been omitted). 
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